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ABSTRACT 

Aims. Multi-wavelengths observations of TeV blazars provide important clues into particle acceleration in relativistic jets. 

Methods. Simultaneous observations of Mkn 421 were taken in very high energy y-rays (> 200 GeV, CAT experiment), X-rays (RXTE) and 

optical (KVA). Multi-day RXTE observations are also presented, allowing for detailed modelling of the spectral variability. 

Results. Short timescale (as 30 mn) variations in VHE y-rays are found, correlated with X-rays but not with the optical. The X-ray spectrum 

hardens with flux until the photon indices saturate above a threshold flux « 10~ 9 erg s _I cirr 2 . The fractional variability decreases from X-rays 

to optical as a power-law with F TO oc E° - 24 * - 01 . The full spectral energy distribution is well-fitted by synchrotron self-Compton emission from 

cooling electrons injected with a Maxwellian distribution of characteristic energy y\,. Fluctuations in the injected power with P;„j oc y A b explain 

the observed variability. 

Conclusions. The spectral saturation and the power-law dependence of the fractional variability are novel results which may extend to other 
TeV blazars. The ability of Maxwellian injections to reproduce the observed features suggest second-order Fermi acceleration or magnetic 
reconnection may play the dominant role in particle acceleration. 

Key words. Acceleration of particles - Radiation mechanisms: non-thermal - Galaxies: active - BL Lacertae objects: individual: Mkn 421 - 
Gamma rays: observations - X-rays: galaxies 



1. Introduction 

Mkn 421 was the first Acti ve Galactic Nuclei (AGN) to be es- 
tablished as a TeV emitter JPunch et al!l992|) . Such very high 
energy (VHE) y-ray emission is thought to arise from parti- 
cles accelerated in a relativistic jet directed along our line- 
of-sight. Continuous improvements in Atmospheric Cherenkov 
Telescopes (ACT) have resulted in new detections of sources, 
shedding new light on the underlying processes, notably in 
conjunction with X-ray observatories. Mkn 421, Mkn 501 and 
1ES 1959+650 are the brightest VHE extragalactic objects de- 
tected so far, with VHE flares reachin g levels m any times that 
of the Crab flux (see e.g. Konopelkoetal .120031) . Bright objects 
like these can be studied with ACTs on smaller time scales than 
fainter objects since there are more photons. 

Variability in Mkn 421 is present at all wavelengths and 
correlations during high states have been found from millime- 
ter wavelength observations up to VHE energies suggesting a 
single population of particles may be responsible, if only partly, 
for radiation over most of the spectrum. The flux doubling 
timescales vary widely but are correlated with energy as the 
highest energies vary faster. 

The spectral energy distribution (SED) has the typical 
double-humped blazar feature in vF v representation, where the 
peaks are separated by 8 to 10 orders of magnitude in en- 



ergy. The lower component is attributed to synchrotron radi- 
ation from relativistic electrons/positrons, and peaks at a fre- 
quency vs which ranges from w 0. 1 ke V in a low activity state 
to a few keV in flaring cases , being clearly correlated with the 
intensity dTanihata et alJl2004T) . The higher component can be 
interpreted as inverse Compton scattering of leptons on the syn- 
chrotron photons themselves. An alternative to leptonic models 
are the so-called "hadronic models" where the relativistic jet 
has a significant component of relativistic protons emitting X- 
rays and y-rays via synchrotron, photo-prod uction, or piqn de- 
cay from proton-proton interactions (see e . g. ICostamantel 2004 
for a recent review). 

Mkn 421 was in a historically high state in the first months 
of 2001, triggering a multi-wavelength campaign and spawning 
a wealth of observations with unprecedented detection levels 
in most energy bands. The campaign included amongst oth- 
ers the ground-based CAT Cherenkov telescope (200 GeV - 20 
TeV), the RXTE X-ray observatory (2-30 keV) and the 60 cm 
KVA telescope on La Palma. The RXTE follow-up campaign 
was optimised for the Whipple ACT, but some observations oc- 
cured during nights at the location of the CAT observatory. An 
analysis of the 4 ks segment that happened to exactly overlap 
with CAT and its impact on SED modeling, are discussed here. 
Multi-day spectral analysis of the X-ray data is also presented. 
Observations and analysis are found in §2, the results are col- 
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lected in §3. A model of the SED is developped in §4 and the 
interpretation of the results is discussed in §5 before conclud- 
ing. 

2. Observations and data analysis 

The data taking, reduction and analysis for CAT and RXTE 
is described next. The KVA telescope data shown in this pa- 
per have been presented by Sillanpaa et al.l J200ll) where the 
reduction and analysis is described. 

2.1. CAT 

The CAT (Cherenkov Array at Themis) imaging telescope op- 
erated on the site of the former solar plant Themis in the 
French Pyrenees from 1996 to 2003. The detector collects 
Cherenkov light from the particles in atmospheric air-showers 
with a 17.8 m 2 dish. The light is recorded by a high definition 
camera with a 4.8° field-of-view, comprised of a central region 
of 546 phototudes in a hexagonal matrix spaced out by 0.13° 
and of 54 surrounding tubes in two "g uard rings". A full d e- 
scription of the detector can be found in lBarrau et alJ ll998). 

The observations were carried out on the night of 23-24 
March 2001 when the source was very active in X-rays. The 
observations were made in the so-called wobble mode, i.e. the 
source position S is shifted by 0.29° from the field-of-view cen- 
ter. Showers pointing to this position are considered as ON- 
source data. Showers pointing to the symmetrical position S' 
with respect to the center are considered OFF data and have 
been used to monitor the cosmic -ray background. 

The y-ray image analysis, i.e. the method and the cuts, i s 
based on the standard CAT procedure (Le Boh ec et alj fl*998). 
Individual events are compared to abacus of theoretical average 
images in order to derive the primary direction and the energy 
of events. The extraction of the y-ray signal from the remaining 
background events is base d on a novel statistical method using 
a maximum likelihood fit (Khelifi 2001). This method gives an 
estimation of the number of y-rays, N r> /, and of the source po- 
sition, 5 / given the observational informations of each shower, 
the fitted arrival direction, D, the fitted shower impact param- 
eter on the ground, IP and the pointing elevation of the tele- 
scope, 6. The likelihood function, X, is defined by the product 
of the probability density of each event over all events passing 
cuts in the field of view (Ntot): 

Ntot 

£(D\N 7 j,S f ,IP,0) = Y\ p(D\N yJ ,Sf,IP,6) 

(=i 

The probability density p of each event can be separated as 
a function of the shower type, i.e. a y-ray origin or an hadronic 
origin: 

p = N^/Ntot x Py (D\Sf, IP, 9) + 
(l~N 7 j/N T aT)Xph(D\IP,0) 

The term p y is the point spread function (PSF) of the CAT de- 
tector. It is determined from Monte-Carlo simulations of y-rays 
with a power-law distribution in energy and for different bins 
in elevation angles. The exact value of the power-law index is 



not critical because the PSF extension is dominated by the low- 
est energy events. A differential spectrum c/N/c/E oc E" 2 7 was 
chosen. The term pi t is the probability to reconstruct a shower 
arrival direction D for a hadronic shower. 

As the true arrival directions of hadronic events (cosmic 
rays) are isotropic, pi, can be parameterized from data taken in 
fields of view without any known y-ray source. 

The likelihood X, is maximised numerically according to 
the parameters N y j and Sj. For the analysis of the Mkn 421 
data, the source position S / has been fixed to the AGN posi- 
tion in order to achieve a better y-ray extraction. This statisti- 
cal method has been tested on simulations and on data taken on 
the Crab Nebula. No bias on the reconstructed number of y-ray 
events is detectable, whatever the signal strength. 

The significance of the signal is calculated from ratio of 
the maximum likelihood calculated on the data, £.(N y j) to the 
likelihood calculated by assuming that no signal is present, 
-C(N y — 0). The exact formula of the significance estimation 
is: 

o-= ^2xlog(£j(N y j)/£(N y = 0)) 

It has been checked on simulations that this significance esti- 
mation follows a Gaussian distribution of variance equal to 1 
in absence of any signal. 

The characteristics of the cuts described above and this sta- 
tistical method for the signal extraction have been computed 
with a y-ray spectrum of spectral index of T=2.55 (dN/dE oc 
E~ r ), an integral flux comparable to the one of the Crab Nebula 
and at a telescope elevation of 0=60°. About 41% of the y- 
ray events are kept and the expected significance for 1 hour of 
observation is 5.7 cr/ Vlhour. This method gives better perfor- 
mances th an a method based on a c lassical subtraction "ON- 
OFF" (e.g. lDiannati-Atai et alll99^ . 

The determination of the spectra and the integral fluxe s 
are made with the standard CAT method JPiron etalJbOOlh . 
Briefly, the spectrum extraction is based on forward-folding 
method in which a parametrization of the spectrum shape is 
assumed. And the fluxes are calculated assuming a spectral 
shape derived by the previous method. In both determinations, 
the number of y-ray events used for them is extracted by the 
maximum likelihood method described previously. 

After selecting data based on clear weather conditions and 
stable detector operation, a total of 2.3 hours (corrected from 
dead time) of observations are kept, all taken within a elevation 
angle between 70° and 90°. The number of y-rays extracted by 
the maximum likelihood method are 1362 ± 50 were detected. 
The flux was high enough to have significant data points with a 
10 minute sampling. This lightcurve is represented in panel (a) 
of Figure n The integral flux is derived as a function of time 
assuming a non-evolving spectrum throughout the night (as it 
is justified by the study of the hardness ratio in the section lXTV 
A flux decrease of 63% (significance of 2.3cr) within 30 min is 
followed by a doubling of the flux in 30 min again (significance 
of 3.8cr). These variability timescales had not been detected by 
CAT on Mkn 421 before these observations. 

The VHE time scale variability of Mkn 421 has been ad- 
dressed in past observations. The Whipple observatory's 10 m 
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telescope repor ted a y-ray rate halvi ng time scale of ^ 1 h on 
a flare in 1999 JlVIaraschi et aljfl99^) . but this has to be untan- 
gled from effective area effects before it can be assumed that 
the variability is intrinsic. The fastest doubling time scale re- 
ported by Whipple was 15 min Idaidos et al.l [l996). also ex- 
pressed in rate rather than flux, though it can probably be as- 
sumed that effective areas do not double on this timescale. No 
reference could be found from Whipple where the lightcurve 
is expressed in integrated flux. The HEGRA collaboration ob- 
served the same flaring period but due to poor weather could 
not report on this specific episode. A lightcurve derived in units 
of integrated flux above 1 TeV shows a halving timescale o f 
< 1 h for the night of MJD 51939/40 JAharonian et alJl2003t) . 
The halving timescale reported here is therefore consistent with 
other ACT measurements, but it is the fastest variability seen in 
an effective area independent way. 




10 5 991.90 991.92 991.94 991.96 991.98 992.00 992.02 

Time (MJD-51000) 

Fig. 1. Top: CAT lightcurve on the night of 23 March 2001 ex- 
pressed as integrated flux above 200 Ge V. The flux decreases 
by a factor > 2 within 30 min which is the fastest variabil- 
ity seen on Mkn 421 by CAT. Bottom: RXTE 2-10 keV flux 
lightcurve for observation times exactly overlapping or con- 
tained within a CAT measurement. The variability in the VHE 
band is larger than in the X-rays. 



2.2. RXTE 

The X-ray data set used here is publicly available covering 
an outburst in Mkn 421 between MJD 51986 and MJD 51994 
(the whole ObsID has a set of smaller observations that extend 
until MJD 52000 but were not included here). The source was 
detected with a rate that required the bright background model. 

1 ObsID 60145 which can be obtained from the HEASARC website 



The high rate allows to segment the data into segments of a 
few 100 s from which a spectrum and lightcurve expressed in 
erg crrr 2 s _1 can be derived instead of the usual counting rate. 
Hence, the correlation between flux and spectral indices can be 
studied on the same time scales instead of using rate vs hard- 
ness ratios. This can be particularly relevant if, such as the case 
here, the spectrum changes significantly when the brightness 
changes. 

A small part of the Mkn 421 observation campaign of 
RXTE is simultaneous with CAT, corresponding to observa- 
tions taken between MJD 5099 1 .8 and MJD 50992. 1 . This seg- 
ment is analysed such that the RXTE data segments are made 
to match the CAT time segments. This means that the X-ray 
spectral file from which the flux is derived is at most as long as 
the CAT segments from which a flux, hardness ratio or spec- 
trum was derived, the latter being just one big segment lasting 
the whole observation night of 24 March 2001. 

The reduction scheme for the X-ray data set is the follow- 
ing: 

- The STANDARD2 data are extracted in GOODTIME inter- 
vals constrained to be either a segment at most 400 s long, 
or at most as long as the overlapping length of a VHE ob- 
servation by CAT, or at most as long as the whole CAT data 
set used to derive the VHE spectrum. Due to the variable 
active PCU configuration during these observations, this is 
done for each PCU independently. 

- The PCU dependent PHA files are combined using the 
addspec tool weighted by the counts information delivered 
by fstatistic and then the corresponding response ma- 
trices were combined with addrmf . The bright background 
model was used and only the 3^40 PHA channel range was 
kept in XSPEC v . 12 . 2, or approximately 2-20 keV. 

- These segments are then fitted to a power-law and a broken 
power-law in XSPEC with PCU configuration-dependent 
response matrices generated by the ftool pear sp v . 10.1 
and a fixed column density of Nh = 1.6 X 10 2l) crrr 2 ob- 
tained from PIMMS 2 . The chisquare probability P(x 2 ) for 
both models is estimated. For the time averaged spectrum, 
a double broken power-law is used. 

- This yields the flux and the error (defined here as the 
1-sigma confidence level) on the flux reported in the 
lightcurve, in units of 10~ n erg s _1 crrr 2 in the 2-10 keV 
band, as well as spectral parameters for each of these inter- 
vals, namely the lower photon index Fj_, the break energy 
Ebr and the upper photon index Fh- The broken power-law 
fit yielded better fit results than a single power-law in 97% 
of the fits. If the single power-law is preferred with a higher 
P(X 2 ) then both indices are set to the single power-law 
value. Since the data segments do not have similar lenghts, 
the error on the flux can vary but is correlated with the inte- 
gration time (shorter segments have an higher associated 
flux error). A small amount of segments had integration 
times < 200 s and were discarded from the analysis to keep 
those with a meaningful statistical error. 



2 See http://legacy.gsfc.nasa.gov/Tools/w3pimms.html 
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Fig. 2. Parameters derived from the broken power-law fits to 
the observation segments, as a function of MJD. The lines in 
the upper panel indicate the segment where simultaneous ob- 
servations were possible with CAT. Whereas there is a clear 
indication of modulation and a flux dependency of both Fl and 
Th, no such indication is visible for the break energy. 



The lightcurve derived to match the CAT lightcurve is rep- 
resented in panel (b) of Figure [2 The total lightcurve for the 
analysed proposal 60145 is in panel (a) of Figure|5]along with 
the curve indicating the variation of the two photon indices 
Tl and Th over the same period in panel (b). The integrated 
flux lightcurve has also its advantages over a simple count- 
ing rate lightcurve since the rate doesn't properly indicate how 
spectral changes affect the flux. Using the photon indices in- 
stead of the hardness ratio gives the absolute spectral varia- 
tions. The lightcurve shows tremendous variability, with a fac- 
tor of ^ 10 variation in integrated flux on timescales of 1 d. A 
similar lightcurve expressed i n X-ray counts/secon d from this 
same data set can be found in lFossati et alJ (120041) along with 
the Whipple VHE curve though no spectral X-ray information 
was derived in i t. A fraction of the same RXTE data set was 
also analysed bv lNolan et alJd2004 . 



3. Results 

3.1. VHE and X-ray Variability 

On the night of 2001 March 23, the RXTE/PCA and CAT ob- 
servations overlapped exactly for over 2 hours of observation 
time. During this time the source underwent a major amplitude 
variation recorded by both instruments as seen in Figure^ The 
amplitude change was much more dramatic in the VHE band, 
varying by as much as a factor of 3 whereas the X-ray flux 
varied by 60%, implying a difference of a factor of 5 in flux 
variability difference. 

In order to quantitatively look for correlated variability, the 
VHE flux is shown versus the X-ray emission in Figure[3] The 
correlation factor for this data set is r = 0.86. The overall cor- 
relation plot shows a linear relation between both fluxes of the 
form 0(E > 200GeV) oc (12.8 ± 1.8) X F 2 -iok e v- The amplitude 
of the variability in the simultaneous light curves is not large 
enough show a deviation from a linear relation, but in order to 
compare with other results, where F vhe k relationships are 
quoted, a power-law will be used in what follows. 

In order to check any difference between a rising and de- 
caying part of the emission, the decaying and the rising parts 
of the lightcurve are fitted to a power-law, yielding k = 3.8 ±0.6 
for the decaying and k = 2.7 ± 0.7 for the rising parts. Within 
the error bars no clear difference can be established between the 
two. Combining all the points for this transient episode yields 
k = 2.9 ± 0.6. The Whipple obser vations from this pe riod were 
also correlated with the X-rays. In Nol an et al. 1 120041) the VHE 
band was correlated with the ASM rate, where the VHE/X-ray 
correlation is not conclusive. This is not surprising sin ce the 
ASM is less sensitive than the PCA. lFossati et alJ 12004) show 
correlations between the Whipple y-ray rate and RXTE/PCA 
counts with k = 1.3+0.3 for the whole data set and a- = 2.3+0.3 
for an individual transient i.e. on a shorter timescale. The value 
of k — 2.9 + 0.6 derived for the CAT/RXTE correlation on short 
timescales is compatible with that result. Note however that in 
the case of a changing spectrum, the count rate of the PCA is 
not directly propor tional to the 2- 10 keV flux which is derived 
from a spectral fit. Tani hata et al.l J2004I) find k — 1.7 ± 0.3 over 
a week during the 1998 campaign, using the peak synchrotron 
flux determined from ASCA, EUVE and RXTE; and VHE y- 
ray flu xes in Crab units from HE GRA and Whipple. As pointed 
out by Katarzvhski et al.1 d2005l) on the basis of observations of 
both Mkn 421 and Mkn 501, correlating long-term lightcurves 
may give results that are significantly different from short-term 
transient correlations. 

A hardness ratio HR was derived for the VHE data set, with 
HR = "^\\, where N(> 1.) and N(> 0.35) are the number of 



AT(>0.35) 

y-rays with energy above 1 TeV and 350 GeV respectively. The 
length of the segments used to derive the ratio was enlarged to 
a meaningful statistical significance (panel (a) in Fig.0}. The 
HR does not appear to be variable for this episode since it is 
compatible with a constant function of time {P(x 2 ) = 0.13). 

The RXTE data were analysed within these HR time seg- 
ments, from which photon indices F L and r H , the indices for 
the power-law component respectively below and above the 
break energy E\, T in the broken power-law, were derived. The 



Berrie Giebels et aL: Unveiling the X-ray/TeV engine in Mkn 421 



5 



photon index as a function of time is in panel (b) of Figure 0] 
The variations of T L and F H are not compatible with a constant 
(P(X 2 ) = 0.04 and 2.3 x 10~ 5 respectively) since they are cor- 
related with the flux. 




10 11 12 13 14 15 
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Fig. 3. VHE flux versus the X-ray flux for the lightcurve shown 
in Figure^ The correlation factor r = 0.86 indicates that both 
wavebands are strongly correlated. There is an indication of 
a counter-clockwise pattern in this diagram where the VHE 
flux is lower in the rising part than in the decaying phase of 
the episode seen here. The decaying and the rising parts of the 
lightcurve are fitted to a power-law (dashed lines). 



3.2. X-ray and optical variability 

The optical data fro m the KVA telescope presented in 
JSillanpaa et aljEoOlh had a few overlapping segments with 
the RXTE observations. The lightcurves of Fig. [5] show that 
there are no counterparts to the rapid and large variations seen 
in the X-rays, as has been mentioned in other similar multi- 
wavelenght studies. There is hence no evidence for simultane- 
ous correlation between those wavebands, in agreement with 
othe r multi-wavelength camp aigns carried out on this object 
(e.g. lBlazeiowski et alJ l2005). The absence of correlation be- 
tween simultaneous optical/VHE lightcurves is patent i n othe r 
VHE emitting B L Lacs such as Mkn 501 JPetrv et alj EoOO). 
1ES 1959+650 jKrawczvnski et alJ l200l or PKS 2155-304 
jAharonian et all2005l) 

To further investigate the energy-dependent variability, the 
fractional root mean sq uare (rms) variability amplitude F val 
(see e.g. IVaughan et alJl2003l for F var and a derivation of the 
uncertainty o> var in equation B2) is estimated for the X-ray 
and the optical lightcurves. F var is defined by 
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Fig. 4. (a): Hardness ratio from the CAT data set, which is com- 
patible with a constant, (b): Photon indices Tl and Fh on the 
same time scale as CAT where variability is evident for both 
indices due to their correlation with the flux. No correlation is 
exists between the VHE HR and the X-ray indices. 

where s 1 is the variance of the considered time series, cr err is 
the mean square error of the flux measurements and x the mean 
of the flux. 

Figure[6]shows F var as a function of energy along with the 
optical KVA measurement of F val = 7.55 ± 0.71%, showing 
evidence for a power-law behaviour of F vm over four decades 
in energy with F var oc ^O- 24 * 01 The increase in fractional vari- 
abili ty with energ y between 0. 1-10 keV was previously noticed 
by iFossati et alJ j20o'o3) using BeppoSAX data. Their X-ray 
fractional rms is also oc E 1 ^ 4 . The historical EUVE data point 
shown in Fig. [6] taken from the same intensive 1998 campaign 
on Mkn 421, is in very good agreement with the power-law 
trend. 

3.3. X-ray spectral patterns 

It is commonly reported in blazar studies that the X-ray flux 
(or rate) correlates with the spectral index (or hardness ratio), 
indicating that the spectrum becomes harder when the flux in- 
creases. Since the X-ray data analysed here were fitted to a 
broken power law the flux dependency of 3 spectral parameters 
can be checked, namely the photon indices and the break en- 
ergy £bi . Surprisingly the break energy varies very little with 
time and also does not correlate with the flux, a fit to a con- 
stant yields ax 2 of 1F6 for 9 degrees of freedom, or a 23% 
chisquare probability. The distribution of E\, r has an average 
value of Ehi- = 5.9 ± 1.1 keV independent of the flux. 

The photon indices Tl and Th however show some inter- 
esting features where the flux dependency is apparent up to a 
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Fig. 5. X-ray flux (crosses) and KVA optical flux (points), for 
3 nights where observations overlapped. Note that the Y axis 
numbering is valid for both wavebands. The correlation factor 
r = 0.06 indicates that both wavebands are not correlated, as is 
apparent in the lightcurves. When the 2-10 keV flux varies by 
more than 50% no such behaviour is visible in the optical. 



flux F„ 



10 9 erg cm 2 s 1 



after which the indices appear to 



plateau (Fig.0. Menti ons of a possib le hardn ess ratio satura- 
tion were made in lGeorge et alJ 0988I) and lNolan et alJ y004) 
but the datasets were insufficient to establish this behaviour. 

A simi l ar an alysis was performed on Mkn 501 by 
IXue & Cull J2005I) who derived the same broken power-law 
spectral parameters for that source and also found that the two 
photon indices show a spectral hardening and then level off at 



a flux of F„ 



5 x 10 10 erg cm 2 s '. 



The authors used a func- 



tion consisting of a constant plus an exponential to describe the 
trend. This functional form does also fit quite well the Mkn 421 
data shown here, as seen in the dashed lines in Figure^! 

The two photon indices are correlated (r = 0.99) and a 
linear fit to the data set yields Th °c (0.96 + 0.02) x Tl showing 
that the dependency is very linear, unlike the flux dependency 
found previously. The ratio between the two indices remains 
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Fig. 6. Fractional rms as a function of observed energy, for the 
RXTE and KVA data from this campaign. Also a EUVE mea- 
surement is sh own as an open circ le, from a MWL campaign 
on Mkn 421 frakahashi ~t all2000h in 1998. The data are well 
described by a power law with a slope oc E 1 ^ 4 . 



therefore constant throughout the whole observation sequence 
here, and is independent of the measured intensity. The part of 
the X-ray spectrum here appears to be pivoting around the 6 
keV break energy, which is probably more an indication of a 
spectral curvature than a feature in the spectrum. 

3.4. Resolved flares 

As noted by |2004) the shortest rise or decay times in 
Mkn 421 can be as short as 1000 s, but clearly resolved ris- 
ing, peaking and decaying structures are not easy to define in 
the lightcurve due to the erratic exposure. There are however 
two rising and decaying structures visible in the intervals MJD 
51987.0-51987.7 (flare A) and MJD 51992.6-51993.0 (flare B) 
with different amplitudes. Flare A's flux is < F p whereas Flare 
B is > F p . This provides an interesting test for the spectral in- 
dex behaviour: both flares start from a different baseline level, 
double their flux, and don't exhibit the same spectral evolution: 
only the fainter flare shows a clear spectral evolution (Figure 
|HJl. The fact that the spectral index varies more when the flux is 
below the F p limit, both in the overall lightcurve as well as in 
resolved flares, will be discussed further in 34.41 

3.5. Time averaged spectrum 

The CAT spectrum is derived from data between 300 GeV 
and 5 TeV, taken at zenith angles taken in the range 0° - 
20°. The most probable spectrum for this data is 4j? = 

<D x £-»-/sxiog 10 £/iTev with o„=(13 + 0.85.„ fl , ± 2.76 JVJ/ ) x 

lO^'cnrV'TeV- 1 , a=2.81 ± 0.09.„ fl , ± 0.07 sysl and /3=0.79 ± 
0.27, (fl/ ± 0.03 !vsr . The curvature term is significant at 2.7cr. 
The same data set can be fitted to a power law with an ex- 
ponential cutoff yielding a similar likelihood with a cutoff en- 
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Fig. 8. Spectral variations for a resolved flare where the flux is < F p (left panel) and > F p (right panel). In the former case, both 
indices harden from their baseline value, peak later than the flux, and soften again with the characteristic hysteresis loop. The 
spectral indices are highest well after the flux peaks. In the latter case, no significant spectral variability is found when the flux 
doubles. The lines represent a fit of the flares to a lognormal function above a baseline modelled by an exponential decay (dashed 
line). These functions help to locate the extrema of the lightcurves but are not used to derive any physical parameters. 



ergy of Eq — 1 .4 + 0.5 TeV. The spectral cur vature here is 
consistent with HEGRA jAharonian et alJl2003t) on Mkn 421 
during these outbursts (though with a different time coverage) 
where a spectral cutoff at Eq = 3.4 ± 0.5 TeV at the 2.4<x level 
was reported. Th e Whipple data also show a curved spectrum 
( Kre nnrich et alJl20 02) compatible with an exponential cutoff 
at Eq = 4.3 + 0.3 TeV as well as the fact that there is no ev- 
idence for variation in the cutoff energy with flux. The differ- 
ence in the value for Eq from these experiments can be due to 
systematic effects (the errors quoted here are statistical only, all 
experiments give a systematic error on Eq of ~ 20%) and to the 
different observation times. 

In order to derive a time-averaged X-ray spectrum consis- 
tent with the VHE spectrum, the RXTE data are reduced to ob- 
tain PCU-dependent PHA files constrained to start at the begin- 
ning of the first CAT run, and to end at the last run. The spec- 
trum is fitted to a double broken power-law in XSPEC which 
yields ax 2 of 52.1 for 50 degrees of freedom, with an excellent 
chi-square probability of P(x 2 ) = 0.39. The spectral parame- 



ters are T L = 2.20 ± 0.06, T M = 2.37 ± 0.10, T H = 2.64 ± 0.04, 
Ebr.i = 6.2 + 0.1 keV and £br,2 = 13 ± 1 keV where Tm is the 
photon index between the two break energies £br,i and £br,2- 
Hence, there is significant curvature in the time-averaged X- 
ray spectrum, which is in agreement with earli er observations 
of Mkn 421 by tel escopes such as ASCA dTanihata et alJ2 004) 
and BeppoSAX dFossati et alJ l2000bt iMassaro et alJ l2004 . 
The 2-10 keV flux is * 1.34 x 10~ 9 erg cm~ 2 s _1 . 

The main facts that will be discussed now can be summa- 
rized as follows: 



- The X-ray and VHE y-ray fluxes are highly correlated 
within a time window of 3h. 

- The optical flux rms, not its mean value, is correlated to the 
higher frequencies through its fractional RMS. 

- The X-ray spectrum hardens with increasing flux only up 
to some characteristic level F p after which it levels off. 
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Fig. 7. Dependence of the three spectral parameters on the flux, 
all derived from fits to a broken power law. The indices are seen 
to be levelling at some characteristic flux, whereas the break 
energy does not appear correlated with the flux. 



4. Synchrotron self-Compton modelling of 
Mkn 421 

VHE emission from BL Lac objects is typically interpreted as 
due to radiation from particles accelerated to high energies in a 
relativistic jet. The nature of the particles, their location(s), the 
acceleration mechanism(s), the radiative processes involved, 
etc, are all subject to debate. The simplest view is that of syn- 
chrotron self-Compton emission from electrons (positrons) lo- 
cated in a single, homogeneous zone. This view is adopted in 
the present work and confronted to the results accumulated in 
the previous section. 

4.1. Constraints from the variability 

Several simple constraints on the physical conditions in the 
emission zone may b e obtained from the observations (e.g. 
iGhisellini et al.lll996l) . The observed variability timescale sets 
an upper limit to the emitting zone. The VHE doubling 
timescale of 30 minutes (Fig.^ implies 

R6~ l < 5.4 x 10 13 cm (2) 



where R is the characteristic size of the region (assumed spher- 
ical in what follows) and 6 is the boost factor from the bulk 
relativistic motion of the jet. The correlated X-ray and VHE 
emissions indicate that the emitting particles are likely to be co- 
spatial, and hence that the emitting volumes are similar. This 
is consistent with the hypothesis of a homogeneous emission 
zone. 

The detection of y-rays in the 0.3 - lOTeV range implies 
that the opacity to pair production r yy within the emission re- 
gion cannot significantly exc eed 1 , independently of the emis- 
sion mechanisms. Following iTavecchio et alJ i 19981) the opac- 
ity Tyy iS 

_ o-T l_dj 1 / 6 2 \ 

Tyy ~ 5 he R 6\l+z) { '(l+z) 2 }- ( ) 

In the above, T yy is the opacity seen by VHE photons with 
an observed energy E y . The VHE photons see a flux den- 
sity F(E r ) of photons above the threshold for pair production 
E, = (m e c 2 ) 2 /E y , which occurs in the UV to X-ray range (ob- 
server frame). If the emission is co-spatial, the measured X-ray 
data can be used to estimate the flux density that VHE y-rays 
may interact with. Imposing that r yy < 1 for photons with ob- 
served energy E y - lOTeV and using Eq. [2 to constrain R 
yields 

6 > 16. (4) 

The mean energy E of synchrotron radiation from an elec- 
tron with Lorentz factor y is given by E « -£- 2ihy ' <lB where 
B is the magnetic field intensity in the emission zone. If the 
observed X-rays are synchrotron radiation from electrons then 
the maximum observed energy £ mal = 20 keV, corresponding 
to the highest energy bin in the RXTE spectrum, yields 

y max > 1.4xl0 6 ZT 1/2 <T 1/2 (5) 

In the SSC scenario, the same electrons should have enough 
energy to Comptonize ambient synchrotron photons up to the 
VHE regime, requiring 

6 , 

z 7max'W e C Z * lOTeV. (6) 

1 +z 

Combining Eq.|3]and Eq.[6]yields a constrain on the magnetic 
field 

BST 1 <5.2xlO" 3 G (7) 

which, for 6 = 16, would imply B < 0.08 G and y max ^ 3.5 X 
10 s . This is of the order o f the Lorentz factor der ived in a sim- 
ilar way for Mkn 421 by Takahashi et al. ( 1996) who derived 
Tmax > 5 x 10 5 but lo wer than y max » 10 7 found in a similar 
way in 1ES 1959+650 dGiebels et alJl2002l) . Accommodating 
a higher magnetic field would require large relativistic speeds 
(boost factor 5). 

The relative strength of the inverse Compton (IC) contribu- 
tion relative to the synchrotron 'bump' depends on the particle 
and magnetic field densities. Higher densities increase the IC 
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emission while a higher B increases the synchrotron contribu- 
tion. In the Thomson regime, the relative luminosities are given 
by 

Lc _ (vF v ) s (d L /R) 2 /(c6 4 ) 



> 810 6 



16 



(8) 



5 2 /8tT 

where B and R are replaced by their maximum allowed value 
(Eqs.[5]and[7J, (vF Y )s ~ 2 x 10~ 9 erg s _1 cm -2 is the observed 
synchrotron spectral flux, di w 4 10 26 cm is the luminosity dis- 
tance, and the numerical application takes 6 — 16 (Eq.0J. The 
small radius (high particle density) and low magnetic field im- 
plied by the constraints lead to a towering IC component, while 
both components are observed at roughly the same spectral lu- 
minosities. Therefore, this leads to values of 5 closer to m 40 to 
match the relative luminosities. 

In Mkn 421, inverse Compton emission at VHE energies 
occurs in the Klein-Nishina (KN) regime, where the electron 
gives away nearly all its energy to the photon. This regime is 
reached when the boosting electon has a Lorentz factor y such 
that ye' s > 1 where e' s is the energy of the target photon (in 
units of m e c 2 ) in the jet flow frame. With synchrotron peak- 
ing at w 500 eV then e' s « 10 _3 5 _1 and IC on electrons with 
y > 10 3 <5 occurs in the KN regime. The IC cross-section then 
decreases dramatically so that the luminosity in Eq. [H]is over- 
estimated. Nevertheless, the conclusion that a high value of 6 
is needed to fit the IC component should still hold: in the KN 
regime, the e lectrons will lead to a jet-fr ame peak IC luminosity 
e' lc « l/e' s llGeorganopoulos & Kazanasll2003h . which trans- 
lates into 5 > 30 since e' IC as 10 6 <5 _I (corresponding to about 
0.5 TeV in the observer frame). 

Hence, the VHE variability timescale and the SED imply 
large boosts 6 Z 30, R ~ 10 15 cm and B ~ 0.1 G. 

4.2. The standard SSC model 

The validity of the above limit on 6 is investigated using a 
numerical code calculating the homogeneous SSC emission 
from a given distribution of electrons in spherical geome - 
try tfoneset all 1 19741 lGouldlll979t iBand & Grindlavlll985l) . 
Radiative transfer is treated as in iKataokaetalJ i 1 9991) . The 
Bessel functions, need ed for the synchrotron em ission and 
absorption coefficients jRvbicki & Lightm an||l979 |). are com- 
puted using the Chebyshev expansion of MacLeod ( 1996). The 
Compton emissivity is computed using the Ijonesl (11968) ap- 
proximating kernel, which takes into account the reduction of 
the cross-section in KN regime and assumes an isotropic pho- 
ton field. Absorption of VHE y-rays due to pair production with 
the infrared extragalactic background light (EBL) during their 
propagation is take n into account. The pai r production cross- 
section is given in [Gould < fe Schr pderl lll967h . The distribution 
and normalisation of the EBL are set to the upper limits found 
by the H.E.S.S. collaboration using observations of more dis- 
tant TeV blazars jAharonian et all2006l) . 

The shape of the electron distribution is assumed to be a 
power-law with an exponential cutoff dN y oc y~ s exp(—y/yt,)dy 
from y m ; n = 1 to oo. The observed SED (Fig. [9J is power- 
law like from radio to optical with a spectral index or ~ 0.35 
(S v oc v~ a ). If the emission is synchrotron then the underlying 



electron power-law has an index s » 1.7. The decrease in spec- 
tral luminosity at X-rays implies a peak (in vF y) synchrotron 
frequ ency somewhere in the extreme UV range dTanihata et alJ 
2004J). The corresponding y^ is univocally determined by ad- 
justing to the data once B and 6 are set. It is typically of order 
10 5 as y b < y max * 5 lO 5 ^ 1 / 2 ^ 2 (Eq.0with B = 0.2 G 
and 5 = 40). Besides the shape of the particle distribution, the 
model parameters are the magnetic field B, the radius R, the 
boost 5 and the particle density n. The first two are set to their 
maximal value allowed by constraints, depending on 5 (Eqs.|2] 
and[7Ji. This ensures IC emission is maximised. The only free 
parameters left are 5 and n. Fits are attempted starting with 
5-16 and going upwards. Adequate fitting of the SED does 
require values above 6 ~ 40 as argued in 34. II 

The information gained by such modelling is limited. The 
particle distribution adjusted to fit the observations has s*1.7 
when values of 2-2.2 are the c anonical values preferred theory- 
wise ( Achterberg et al. 2001). Using s ~ 2 would require 
Tmin ~ 10 3 {i.e. a limited range of non-thermal particles) in 
order not to overestimate the IR spectrum and enhance the HE 
y-ray spectrum above the EGRET measurement. A high y^ is 
not improbable as y m j n ~ F 2 is expected for relativistic shock 
acceleration. Yet, one parameter is frozen (the index s) only to 
thaw another one (y m j n ) at the price of reproducing less of the 
SED: the exercise is not particularly enlightning. 

The largest limitation is that such simple modelling is, by 
construction, unable to address the wealth of variability infor- 
mation gathered in the previous sections. Injection and cooling 
of particles has to be explicitely taken into account and com- 
pared to the spectral variability. 

4.3. Evolved SSC model: steady-state 

Particles are now continuously injected in a homogeneous 
sphe re, losing energy to radiation until they escape after a tim e 
f esc iMastichia dTs & Kirkll997HChiaberge & Ghisellinill 19991) . 
The evolution of the particle distribution N(y, t) is given by the 
diffusion equation 



^4[». + ««(r, rtl + ^ + fi(r) 

ot dy f esc 



(9) 



where y represents the radiative loss term to synchrotron and 
IC radiation, N/t esc is the energy-independent escape term and 
Q(y) is the injected particle distribution. The radiative losses y 
are computed from the synchrotron and IC power emitted by 
each electron, taking the Klein-Nishina dro p in cross-section 
into ac count. The equation is solved using t he lchang & Cooped 
Jl970h algorithm dPark & Petrosianll 19961) . In all that follows 
the parameters are set to 5 = 40, B = 0.2 G and R = 2x 10 15 cm. 

In the above, cooling will be unimportant if the particles 
escape before emitting a significant fraction of their energy. 
Taking into account only synchrotron radiation, the critical en- 
ergy above which cooling is important is 



Tcool 



3X10 5 SolCcfesc/tf)" 1 - 



(10) 



Low-energy particles radiate away most of their energy if they 
are given the time. It can also be inferred from this equation 
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that light travel time delays must be taken into account when 
investigating the emission of electrons with y X 3 x 10 5 , as 
these should cool on a time £ R/c according to Eq. 1 101 above. 

If the escape time is of order R/c then the steady-state dis- 
tribution is essentially the injected one because none of the par- 
ticles cool significantly before escaping (y coo i ~ y max )- The 
situation is identical to 34.21 However, the gyroradius Rl * 
10 10 76fio2 cm °f TeV electrons is very small compared to R 
(Eq. |3, so it is more likely to have the electrons diffuse for 
some time in the emission region rather than see them escape 
on a direct flight. If f esc » R/c the impact of cooling on the 
SED will be important. 

Classically, with a power-law injection of particles Q oc y~ s 
and j cc y 1 (synchrotron and Thomson IC), the steady-state 
(f — > oo) distribution is a broken power-law with an index 
s + 1 above y coo \. If, as diffusive shock acceleration suggests, 
s — 2 up to some y max , the cooled spectrum has 5 = 3 and a 
flat SED from hv coo i « 13 d^B^ict^/R) ~ 2 keV up to /zv max 
(>20 keV, ^4~TV A flat SEP hardly fits the optical to X-ray 
range: Tanihat a et alJ (2004) find the emission is not flat but 
peaks in the UV. In addition, the emission below v coo i would 
still be a power-law with a = 0.5, inconsistent with observa- 
tions. 

The core problem in reproducing the SED with cooling 
power-laws is the large range of injected energies. There are 
too many low-energy electrons producing too much flux at op- 
tical frequencies and below. One remedy is to arbitrarily harden 
the injection to s = 1 .7 as done in 34.21 fine-tuning jf, so that 
the optical to X-ray transition is reproduced. As previously ar- 
gued, this is not very satisfying. Another is to keep s — 2 but 
increase y m ; n so that high-energy particles carry more of the to- 
tal injected energy, cutting off the low frequency emission. In 
the strong cooling regime where y m \ n > y coo i there is no un- 
cooled component in the distribution. 

Injected particles over a limited energy range was found 
to work very well. In the limit of a mono-energetic injection, 
the steady-state particle distribution yields a power-law of in- 
dex 1 .7 for the chosen parameters, in excellent agreement with 
the observed SED. For synchrotron and Thomson cooling, the 
steady-state is well-known to be a power-law with an index 2. 
Here, inverse Compton cooling occurs in the less efficient KN 
regime, leading to a hardened equilibrium solution. That the 
mono-energetic injection naturally yields the right index with- 
out fine-tuning makes this solution extremely attractive. 

In the following, relativistic Maxwe llian injections 
JSauge & Henril2004tlKatarzvriski et all2006l) 

Q(y) = Ky 2 exp(-2y/y b ) (11) 

were adopted instead of strictly mono-energetic injections or 
power-laws over very limited ranges 3 . The results are identical 
since the Maxwellian is strongly peaked at its characteristic en- 
ergy jb. Maxwellians may arise from second order Fermi scat- 
tering on plasma waves or from acceleration during magnetic 
field reconnection (see §5). 

3 Note that this is not an s = 2 accelerated power-law with expo- 
nential cutoff, y~ 2 exp(-y/yi), as previously discussed. 



Figure [5] shows the result of a fit to the simultaneous opti- 
cal (KVA), X-ray (RXTE) and VHE (CAT) data, along with 
NED archival data and the time-average d spectral mea sure- 
ment from the third EGRET catalog IHartman et alJl9 99). The 
contribut ion of the host galaxy a t optical/IR was modelled fol- 
lowing Katarzvrisk i et alJJ20Q3h as a blackbody of temperature 
r=3500 K, normalised to 7xl0~ 25 erg s~' cirT 2 at its peak. 
The injected Maxwellian distribution has 7^=1.6 x 10 5 with 
K — 3.3 x 10~ 21 in Eq. ^2 This is equivalent to an injection 
rate of n; n j = 9.7 x 10~ 7 particles crrT 3 s _I or a total injected 
power Pj n j = 4 x 10 39 erg s . 

In the strong cooling regime, the escape time sets the min- 
imum energy up to which particles cool, given by y coo i (tak- 
ing both IC and synchrotron into account). With f esc = 50 R/c 
the synchrotron emission cuts off below optical and is un- 
changed above. The VHE emission is also unchanged but the 
GeV flux is lower as there are fewer low-energy seed photons. 
The adopted escape time is f esc = 5000 R/c so that particles 
cool down to y ~ 10 before escaping, emitting low-frequency 
radio-waves. 

This radio emission is self-absorbed in a one-zone model. 
If the radio emission zone is more extended (keeping the same 
B), the density of particles may be low enough to forego self- 
absorption. In general, the interpretation of the radio emis- 
sion from blazars in SSC models requires the addition of an 
addi tional large-sca le component representing the VLBI jet 
JCharlotetalJ2004. 

The electron density in steady-state is ni n jf eS c ~ 325 par- 
ticles cirT 3 . The energy density in particles is about 5xl0~ 2 
erg cm" 3 compared to a magnetic field energy density of 
1.6xl0~ 3 er g cm" 3 . A sub-equipa rtition magnetic field was 
also found bv lMaraschi et alJ Jl999h and appears generic of ho- 
mogeneous models. The total luminosity is « 1.5xl0 46 ergs" 1 . 

In conclusion, a steady-state homogeneous model with a 
Maxwellian-like injection of particles provides a simple de- 
scription of the overall SED, removing the need to set-up a 
power-law with an arbitrary hard index. In the following, this 
picture is confronted with the constraints on variability derived 
in §3. 

4.4. Evolved SSC model: variability 

The observed variability decreases very rapidly with wave- 
length to only a few percent rms at optical (Fig. [5}. X-rays 
should therefore provide the best constraints on the neces- 
sary modifications to explain the variability. The spectra cor- 
responding to the highest and lowest X-ray fluxes measured 
in the campaign are shown in Fig. ^0] This figure shows that 
the simultaneous detection of Mkn 421 by RXTE and CAT de- 
scribed earlier (Fig.|5Jl occurred while the X-ray flux was quite 
close to the maximum observed in this campaign. 

A very good representation of the iow' and 'high' SEDs is 
obtained by varying only the energy of the injected Maxwellian 
by about 50% from y b = 7.6 x 10 4 (low) to 11.4 x 10 4 (high), 
assuming that the source has time to reach steady-state (this 
is relaxed in §4.5). The coefficient K (Eq.ll 1> is kept constant 
to its value of 3.3 X 10" 21 found from fitting the simultaneous 
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Fig. 9. Simultaneous SED of Mkn 421 with archival EGRET 
(grey butterfly) and NED (grey points) data. The simultaneous 
VHE (CAT), X-ray (RXTE) and optical (KVA) data (§2) are 
highlighted in black. Most of the variability in the synchroton 
component occurs at X-ray energies (see Fig.|6jl so that archival 
data below 1 eV, even if non-simultaneous, still provides valid 
constraints. The fit (solid line) is a steady-state homogeneous 
model assuming continuous injection of a Maxwellian distribu- 
tion of electrons with y-53 x 10 4 at a rate n lx ^ =9.7x 1CT 7 par- 
ticles cirr 3 s _1 and cooling via synchrotron and self-Compton 
emission. The contribution of the host galaxy to the optical- 
IR (dotted line) is modelled using a blackbody spectrum as in 
iKatarzvhski et alJ d2003l) . Values for the other parameters are 
6 = 40, B = 0.2 G, R = 2 x 10 15 cm and f esc = 5000 R/c. The 
full line takes into account pair production on the EBL at VHE 
energies and synchrotron self-absorption at radio frequencies. 
The dashed line shows the unabsorbed model. Cooling in KN 
regime of the injected electrons leads to a steady-state distribu- 
tion with the right slope to reproduce the radio to optical slope. 

data. The injected power during the RXTE observations varies 
from I.lxl0 39 -5.7xl0 39 erg s _1 . Hence, the injection behaves 
'thermally' with n; n j oc -y 3 or fi n j oc y A b . The characteristic jb is 
typically set by the balance of the radiative timescale with the 
acceleration timescale f acc but it is unclear how this relates to 
fi n j without a detailed model for the acceleration process. 

Varying Pj n j or n; n j, keeping yi, constant, alters the normal- 
isation of the synchrotron and IC emission, but does not repro- 
duce the spectral softening at X-rays with lower flux. Varying 
both jb and n;„j (or fi n j) is necessary to account for the SEDs. 
Adjusting both independently, there being no strong a priori 
reason why n; n j oc y 3 , still leads to this steep dependence. 

No other single parameter change can reproduce the three 
SEDs. It is difficult to fathom why f esc , R, B or 6 would change 
on a timescale R/c in the context of a one-zone model. In any 
case, varying these one at a time does not reproduce the evolu- 
tion. Combined variations cannot be excluded, but the simplest 
solution is that the injected power fluctuates, giving rise to par- 
ticles with characteristic energy jb oc ■ 

The high and low model SEDs bracket the EGRET long- 
term average and have the same spectral slope. The VHE flux 
varies by an order-of-magnitude between low and high SEDs. 



The low VHE flux is consistent with the lowest measure- 
ments reported by Wh ipple and HEGRA jMacomb et al J 19951 
lAharonian et all2002l) . 

Figure ^2 shows the fractional variation / = (Fhigh - 
^iow)/(^high + ^low) between the low and high SEDs. The 
strongest variability (/ ^ 0.5) occurs between 1-100 keV and 
above about 100 MeV (EGRET band). The drop from 0.1 to 
1 MeV is due to the increasing dominance of the IC compo- 
nent. Interestingly, / increases as a power-law from IR to hard 
X-rays, just as observed in the rms variability vs energy plot 
(Fig. [6} but with a shallower slope. In particular, there is no 
feature associated with the synchrotron peak, which is around 
a keV an d moves slightly to l onger wavelengths with decreas- 
ing flux jTanihata et alJ2004 . 

The variability in the IC component (^ 1 MeV) is less 
wavelength-dependent than the synchrotron component, as it 
mirrors primarily the changes in UV to X-ray seed photons. 
With / w 0.3 the optical variations are overestimated, although 
contamination from the host galaxy could reduce the measured 
variability from the expected value. Indeed, such a contribu- 
tion is included in the plotted SEDs. It would produce a fea- 
ture in the rms vs energy plot (~Fig. ll It . which seems unlikely 
since the extrapolation of the observed X-ray rms is spot on 
the optical. However, note that / is not an rms. In principle, 
one should compute a theoretical lightcurve by varying ac- 
cording to some prescription and compare the rms vs energy 
theoretical plot (rather than /) with the observed one. Hence, 
detailed modelling of the rms plot can constrain how Pj n j varies. 
The initial steps are described below but this is left for future 
work. 

In conclusion, changes between high and low X-ray fluxes 
can be modelled as variations in the injected power of the pile- 
up distribution. 

4.5. Towards the variability engine 

Until now, particles have been considered to have time to reach 
steady-state, thereby neglecting any memory the system might 
have of its past condition. More realistically, variations in the 
injections will occur on timescales ^ R/c in order to repro- 
duce the observed doubling timescale. Fluctuations in the en- 
gine power are mirrored in the emission of the VHE particles 
and then propagate down in wavelength as the particles cool. 
The system may not have time to adjust to the new conditions 
before another variation occurs. The SED received by the ob- 
server will depend on the history of fluctuations in injection. 

It is beyond the scope of the present study to consider 
the resulting lightcurves for various prescriptions on how fVj 
varies. Nevertheless, it is possible to illustrate the expected 
change in the limiting case of a response to a simultaneous, 
step-like change in injection throughout the emission region. 
Starting with the steady-state distribution corresponding to 
the high-state (resp. low-state) of Fig. ^3 the injection jb is 
abruptly set to its low-state (resp. high-state) value after some 
arbitrary time. The resulting lightcurves at various energies are 
shown in Fig. [21 While the change in injection is simultane- 
ous throughout the emission zone, an observer would see the 




Fig. 10. Non-simultaneous spectral energy distribution (SED) 
of Mkn 42 1 with the minimum and maximum X-ray fluxes de- 
tected in this campaign shown. Dashed lines are unabsorbed 
spectra. Dotted line is the host galaxy contribution. The iow' 
and 'high' SEDs were fitted by varying the injected power P™ 
in the steady-state homogeneous model, keeping all other pa- 
rameters identical to those used for Fig. [5] (see 34.41 1. The in- 
jected particle energy behaves as jb °c ■ These changes re- 
sult in wavelength-dependent variations between the two SEDs 
(see Fig. II li . At y-ray energies, the SEDs bracket the long-term 
average measurement of EGRET. Variations in VHE flux by a 
factor ~ 10 between low and high states are expected. 

impact of the variation first on the closest part of the sphere 
and see a change on the opp osite side after a time R/dc (e.g. 
IChiaberee & Ghisellinf l999). This light propagation effect is 
take n into account w hen computing the lightcurves, following 
iKataokaetallfeOOOl) . 

As expected, the response is quickest at X-rays and VHE, 
corresponding to electrons with the highest energies. This ex- 
plains the correlations between these bands. The impact of 
the change in injection only propagates down to the optical 
V band after a delay of ~ a day, as previously discussed in 
iMastichiadis & Kirld Jl997t) . This large delay would explain the 
lack of correlations between X-rays and optical bands. The rate 
of change is a few ks for the X-rays whereas it takes a few tens 
of ks for the optical to drop. Sputters in the variability engine 
on a timescale ;$ 10 ks will be smoothed out in the optical 
compared to X-rays, further reducing the amplitude of the vari- 
ations at optical compared to the estimate made in Fig. ^2 At 
GeV energies, the variation is also delayed and very gradual as 
the IC emission results from an integration over the seed spec- 
trum. Variations at GeV could be even smoother in time than 
at optical wavelengths. Quantitative estimates of these aspects 
require a full simulated lightcurve. 

The correlation between the X-ray and VHE lightcurves 
is shown in Fig. ^] for step changes. The finite light travel 
time and slightly different cooling combine to produce an 
asymmetric loop. At the other extreme, a change in particle 
injection could occur at a slow enough pace that the parti- 
cles are in steady-state, at least at those energies. Looking at 
Fig. this will typically require variations in the injection 



Fig. 11. Fractional flux variation between the low and high 
model SEDs in Fig. [K)] The dashed line shows the expected 
dilution in IR-optical variability if the host galaxy contribution 
is taken into account. The dotted line represents the energy de- 
pendence of the rms variability found in Fig. [6] The fractional 
variation of the model (which is not an rms) has a shallower 
slope. The peak variations do not occur at the synchrotron peak 
but at the synchrotron frequency corresponding to the injec- 
tion ji,. IC becomes dominant above 1 MeV, producing a drop 
in variability. The IC variability is less wavelength-dependent 
since it mirrors changes in seed photons emitted from UV to 
X-ray wavelengths. 

on timescales of many ks. In this case, P m ^ is slowly varied 
between its two extreme values. The flux evolution is plotted 
in Fig. Differences between low-to-high and high-to-low 
disappear as the SED is in steady-state at all times. The corre- 
lation is well approximated by Fvhe « F^ 3 , as expected for 
SSC in the Klein-Nishina regime where Lie is proportional to 
the nu mber of electrons rather than its square (G hisellini et alJ 
1996). This is very close to the long-term trends discussed in 
33.11 Steeper dependences of the VHE flux with X-ray flux 
on short timescales, as found in the CAT data, are not ex- 
plained to satisfaction here. The limited dynamic range may 
bias the dependences found in datasets with a short time base 
( Kat arzvhski et al.ll2005l) . Inversely, these may be able to shed 
more light on the workings of the engine than long-term av- 
erages. As with the energy dependence of the rms, computing 
lightcurves given some prescription for variations in P- m j would 
be needed to properly investigate the issue. 

Figure O shows the evolution of the X-ray spectral in- 
dices T L (2-6.5 keV) and T H (6.5-10 keV) in the model for the 
step-like and slowly varying cases. There is a flattening above 
~ 10~ 9 erg s" 1 cirT 2 when the variation goes from low-to-high. 
There is no such flattening in the opposite direction. The rea- 
sons are not obvious and are related to the interplay between 
injection, cooling and light propagation. Nevertheless, varia- 
tions in index should occur within this loop, in agreement with 
the observations (Fig.0. 

In conclusion, the fits made in 34.41 are used to study the 
response to variations in the injection from the low to the high- 
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Fig. 12. Lightcurve response to a step change in the injection 
of particles from low-to-high state (top) and from high-to-low 
state (bottom). The injected yi, is abruptly changed at t = 
simultaneously throughout the emission region. Light propa- 
gation delays are taken into account and time is counted in 
the observer frame. The X-ray lightcurve is integrated from 2- 
10 keV, VHE is the integrated flux above 200 GeV, GeV is the 
integrated flux between 100 MeV and 10 GeV, V is the optical 
V band flux. Note the rapid response of the X-ray and VHE 
flux, the ~1 day delay at optical and the gradual response of 
the flux in the EGRET band. Correlations between X-ray/VHE 
can be expected. However, the long delays and slower rates of 
change will hide correlations with the optical or GeV fluxes. 

state values and back. The variations are assumed to be either 
instantaneous (step-like) or slow enough for the high-energy 
electrons to reach steady-state, thereby illustrating two extreme 
cases. Fluctuations in X-ray and VHE y-ray fluxes are found to 
be correlated with an index identical to the one derived from 
long time-base observations. The expected variations in X-ray 
photon index vs flux also agree well with the observations. The 
GeV and optical suffer long delays and respond much more 
slowly to sudden changes, which should decrease their vari- 
ability amplitude and hide correlations with X-rays or VHE y- 
rays. Hence, the essence of the observational results garnered 
in §3 can be reproduced with a simple change in injected non- 
thermal power fi n j, assuming the distribution is Maxwellian 
with a peak energy yt, oc P 1 ^ . 

5. Discussion 

The CAT observations presented here confirm that fast vari- 
ability on timescales % 1 hr is a hallmark of VHE emis- 
sion from Mkn 421. Nevertheless, the dataset remains lim- 
ited compared to X-rays, making the identification of a char- 
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Fig. 13. Correlation between X-ray and VHE fluxes plotted in 
Fig-El The black curve shows the correlation for step changes 
from a low steady-state (bottom left) to a high steady-state (top 
right) and back. Arrows indicate the low-to-high and high-to- 
low direction in which the loop is traced. The correlation is not 
symmetric because of light propagation delays. The grey curve 
illustrates the other extreme case in which variations occur on 
a timescale > 10 ks so that the highest energy particles always 
reach steady-state. In this case, the light travel time and cool- 
ing effects have no influence and the change is symmetric. The 
correlation is then found to follow F vhe k F^ 3 , as expected of 
IC in the Klein-Nishina regime. 
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Fig. 14. Correlation between X-ray photon indices Tl (top 
panel) and Fh (bottom panel) with the 2-10 ke V flux. The black 
solid line corresponds to step changes in the injection, while the 
grey solid line is for an injection varying on a timescale » R/c 
(as in Fig. I13i. The direction of change from low-to-high and 
back are illustrated by arrows. The observed indices and fluxes 
shown in Fig.0(on the same scale) fall within the dark loop. 

acteristic timescale uncertain. Rapid X-ray flaring, down to 
sub-hour timescales, appears to be assocated with high X- 
ray fluxes - which is also when Mkn 421 is detected by 
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Cherenkov telescopes llCuil2Q04 . At lower fl uxes, X-ray vari- 
ability seems to occur on longer timescales (Takahashiet al. 
2000). Observations with long time-bases and at lower fluxes 
are needed to characterize better the variability properties at 
VHE energies and compare them to the X-ray properties. 

As with other BL Lac objects, the limit from internal 
absorption of y-rays implies a high relativistic bulk motion 
in the jet 6 > 16. Even higher values of 6 ~ 40 are re- 
quired to model the simultaneous CAT, RXTE and optical SED 
with a single-zone, homogeneous SSC model. These higher 
values are needed to ensure that the emission zone is large 
enough, hence densities low enough, to avoid over-producing 
y-ray emission by inverse-Compton emission. The emission 
region size R — 2 x 10 15 cm is smaller and the magnetic 
field intensity B = 0.2 G slightly higher than that fo und 
in previous works (Takahashi et al. 1996; Mastichiadis & Kirk 
19971: iMaraschietalJ Il999t IChiaberge & Ghisellinil Il999t 
Kat arzvhski et alJ2 003: Kon opelko et alJ2 003). 

High values for 6 contrast with the kinematics deduced 
from radio VLBI observations, which typically imply Doppler 
factors of a few at larger scales. VLBA measurements of 
Mkn 421 realised shortly after the high VHE state reported 
here did not find any new compon ents at the parsec scal e 
after the VHE high state in 2001 JPiner & Edwards! 1200-4 . 
Solutions to this inconsistency include strong jet decelera- 
tion between the y-ray emi t ting z one and the parsec scale 
jGeo rganopou los & Kazanasl 120031) . o r a fast-moving spin e 
surrounded by a slower- moving sheath JChiaberge et alJ2000l) . 
However, high 6 factors are likely inherent to single-zone SSC 
models and more complex, inhomogeneo us model can be ex- 
pected to reduce the discrepancy (Henri & Sauge 2006). 

The spectral variability of Mkn 421 is investigated under 
the assumption of continuous injection of high e nergy elec- 
trons in the emission region. [Mastichiadis & Kirkl II 19971) and 
IChiaberge & Ghisellinil i 19991) had used injection functions 
Q oc y~ s exp(-y/yz,) with s =s 1.7 set in order to reproduce 
both the synchrotron X-rays and the power-law at lower ener- 
gies. This is a harder value than expected for the power-law 
distributions resulting from diffusive shock acceleration. 

A Maxwellian injection naturally yields a cooling distribu- 
tion with an index s x 1.7, forfeiting the need to arbitrarily 
set an injection index. Pile-ups arise naturally from stochastic 
acceleration (second-order Fermi acceleration) on plasma wave 
turbulence and can be as effective in producing hig h energy par- 
ticles as first-order Fermi acceler ation at shocks ( Schlic keiseJ 
11984 Il985t lAharonian etai] [l986). Quasi-monoenergetic in- 
jections may also be expect ed from reconnect ion in Poynting- 
flux dominated jets (e.g. ISikora et alJ 120051 and references 
therein). 

The substantial X-ray spectral variations are well repro- 
duced by a factor * 5 change in injected power. The charac- 
teristic energy of the particles yi, varies as PJ. . A thermal-like 
correspondence between total injected power and characteris- 
tic energy of the VHE particles is surprising: even if the ac- 
celeration process results in a pile-up distribution, why should 
the injection rate vary as y? ? Whether this is a coincidence 
or something which to be addressed by acceleration theory re- 
mains to be decided. 



The injected particles have a characteristic yi, ~ 10 5 , which 
corresponds to the minimum energy for which electrons cool 
in a time $ R/c, and which also correponds to the maxi- 
mum energy that can be reached in the acceleration process 
before cooling dominates. This suggests that f acc ~ f coo i ~ R/c. 
Symmetric rise an d decay lightcurves are e xpected in this case 
at X-ray energies (Chiaberge & Ghisellini 1999), as observed 
in Mkn 421 JCul2004l) . 

The higher yi, hints at a decreasing f acc , hence a more ef- 
ficient particle acceleration mechanism in the high state. This 
also explains the saturation of the X-ray spectral indices above 
a certain flux. The index varies less when the synchrotron peak 
reaches the spectral band. Such a saturation has also been ob- 
served in Mkn 501. A detailed comparison between these and 
other sources could provide simpler diagnostics of accelera- 
tion at the highest energies than hysteresis loops (e.g. Fig. ll3> . 
Both clockwis e and counter-clockwise loops have been seen 
in Mk n 421 fTakahashfiilDIIIIl iFossati et alJboOOht Icul 
2004). The choice of spectral bands and the base-line emission 
may significantly affect the resulting loops and make them am- 
biguous to interpret as isolated cooling/heating events, whereas 
an overall flux - index plot for the whole dataset may provide a 
broader, simpler picture to address. 

Likewise, the power-law dependence of the fractional vari- 
ability with energy from optical to X-ray energies may offer 
new insights into the workings of the engine. It does not im- 
ply correlations between the lightcurves at all these energies 
(i.e. the phase information can be lost). Indeed, the relation- 
ship between rms variability and energy seems unlinked to the 
flux level. In the present model, the engine must somehow be 
able to wag the high energy tail of the particle distribution in 
the right way so that the variation propagates down in energy 
with a tightly-constrained decreasing amplitude. Theoretical 
lightcurves, obtained by varying P; n j, should be computed and 
their rms compared to the observations. Interestingly, the slow 
response of the GeV flux to step changes in injection sug- 
gests strong variations at X-ray and TeV energies which may 
not be accompanied by similar variations at GeV energies. 
Characterizing the amplitude and correlation with other bands 
by detailed modelling will be useful for GLAST. 

6. Conclusion 

Simultaneous observations carried out with CAT and RXTE 
confirm the correlation between X-ray and VHE emission and 
the fast sub-hour variability at VHE energies. A saturation of 
the X-ray spectral indices with X-ray flux and a power-law re- 
lationship between fractional variability and energy (from op- 
tical to X-rays) are found. These results may extend to other 
TeV blazars, and thereby provide strong, universal constraints 
on models of high energy emission from relativistic jets. 

The SED and variability are well reproduced by a single- 
zone, homogeneous model in which a Maxwellian distribution 
of electrons with jb ~ 10 5 are continuously injected. Such 
an injection would favour second-order Fermi acceleration or 
neutral point acceleration over the more widely studied shock 
acceleration. Cooling by synchrotron self-Compton emission 
in the Klein-Nishina regime naturally leads to a steady-state 
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distribution of electrons with the correct slope to explain the 
emission down to radio frequencies. The X-ray spectral satu- 
ration and power-law rms spectrum are reproduced by a factor 
5 change in the power P; n j oc y l J 4 of the injected Maxwellian. 
The acceleration process is probably responsible for this de- 
pendence, while variations in P; n j depend on the workings of 
the internal engine. Detailed lightcurves, particularly at y-ray 
energies where the information is sparse, can therefore provide 
important insights into the mechanisms at work in blazars. 
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